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Abstract 

In  this  paper,  we  present  the  floating  random-walk  solution  of  a  transverse  magnetic 
electromagnetic  problem.  In  a  previous  work,  we  have  developed  a  floating  random-walk 
algorithm  suited  to  the  solution  of  Maxwell-Helmholtz’s  equations  in  heterogeneous  problem 
domains.  The  task  of  developing  an  approximate  expression  for  the  finite-domain  Green’s 
function  for  Helmhotz’s  equations  in  heterogeneous  problem  domains  had  been  accomplished 
with  the  help  of  a  novel  use  of  iterative  perturbation  theory.  We  had  applied  the  algorithm  to  the 
skin-effect  and  transverse  electric  electromagnetic  problems  in  both  homogeneous  and 
heterogeneous  problem  domains.  In  this  work,  we  extend  this  algorithm  to  a  transverse  magnetic 
electromagnetic  problem  and  present  the  preliminary  results  for  a  homogeneous  benchmark.  We 
believe  that  with  additional  development  this  algorithm  will  be  useful  for  the  electromagnetic 
analysis  of  IC  interconnects  at  high  frequencies. 

Index  Terms-  Floating  random-walk  algorithm,  IC  interconnects,  Green’s  function,  Iterative 
perturbation. 


1.  Introduction 

Recent  advances  in  digital  IC  interconnect  technology  have  led  to  multi-GHz  operating 
frequencies.  At  such  high  frequencies,  electromagnetic  effects  such  as  skin-effect  loss,  frequency- 
dependent  inductance  and  capacitance,  slow-wave  substrate  coupling,  distributed  transmission¬ 
line  propagation,  and  high-frequency  radiation  play  a  very  important  role  in  the  electrical 
performance  of  circuits.  Fortunately,  electromagnetic  effects  in  nature  are  described  by 
Maxwell’s  equations.  As  a  result,  the  semiconductor  device  industry  has  placed  a  considerable 
emphasis  in  developing  numerical  solvers  for  Maxwell’s  equations. 

In  a  previous  work  [1],  we  have  developed  a  floating  random-walk  [2]  algorithm  for  the  solution 
of  Helmholtz  equation  in  heterogeneous  problem  domains.  This  algorithm  requires  no 
discretization  meshing  of  either  the  volume  or  the  surface  of  the  problem  domains  and  hence 
memory  requirements  are  expected  to  be  orders  of  magnitude  less  than  discretization-based 
methods  such  as  the  finite-element  [3]  or  the  method  of  moments  [4].  The  algorithm  is 
completely  parallelizable  and  a  linear  increase  in  computational  speed  is  achieved  with  the 
increase  in  the  number  of  processors.  But,  before  we  describe  the  details  of  this  floating  random- 
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walk  algorithm,  we  will  first  describe  the  transverse  magnetic  (TM)  problem  that  we  are  going  to 
solve  in  this  work. 


2.  The  transverse  magnetic  problem 

For  the  TM  problem  [5],  the  electric  field  and  the  direction  of  propagation  are  in  the  xy  plane, 
while  the  magnetic  field  is  in  the  z  direction.  The  electric  field  components  have  possible  spatial 
variations  in  the  xy  plane;  therefore,  E  =  Ex(x,y)ex  +  Ev(x,y)ey  where  ex  and  ev  are  unit 
vectors  in  the  x  and  y  direction.  The  resulting  field  equation  can  be  written  as 

V2E  +  /2E-V(VE)  =  0,  (1) 

where  the  propagation  constant,  y  is  given  by  the  relation 

y2  =  /Lie or  -  jcoafi.  (2) 

Equation  (1)  is  not  of  the  vector-Helmholtz  type;  the  components  of  the  electric  field  are  coupled 
because  of  the  V(VE)  term.  To  achieve  uncoupling  of  these  components,  we  transform  the 
problem  to  a  vector  potential  formulation.  The  four-vector  potential  (cp,  A)  is  a  complete 
description  of  the  electromagnetic  field,  where  cp  is  the  scalar  potential  and  A  is  the  vector 
potential.  In  the  divergence-free  Coulomb  gauge  (i.e.,  V-A  =  0),  each  component  of  A 
independently  satisfies  its  own  scalar  Helmholtz  equation 


V2A  +  y2A  =  /uJ.  (3) 

Above,  J  =  crV (f>  and o  =  afx,v)  +  jcos{x,y),  where  J,Q),a,£  and  //are  the  current  density, 
frequency,  conductivity,  permittivity  and  permeability  respectively.  As  a  result,  we  need  to  solve 
a  scalar  Helmholtz  equation  for  each  component  of  A.  Harmonic  time  dependence  with  a 
complex  propagation  constant  has  been  assumed  in  deriving  equation  (1)  to  equation  (3).  But 
before  we  describe  the  specifics  of  our  benchmark  problem,  we  will  briefly  describe  our 
previously  developed  iterative  perturbation-based  floating  random-walk  algorithm. 

3.  The  floating  random-walk  algorithm 

The  objective  of  the  floating  random-walk  algorithm  is  to  solve  for  equation  (3),  by  solving  the 
scalar  Helmholtz  equation  for  each  of  the  components  of  the  vector-potential.  The  Green’s 
function  equation  for  equation  (3)  is  given  by 

V/G( r  |  r0)  +  y2(r)G(r  ro)  =  A(r-ro),  (4) 

where  G(r|ro)  is  the  2D  Green’s  function  at  position  coordinate  r  with  a  Dirac  delta  function 
applied  at  rD .  If  we  impose  Dirichlet  boundary  condition  onG(r  |  ro),  we  haveG(r  |  ro)  =  0 ,  for 
r  e  S,  where  S  is  a  circle  of  radius  R.  If  we  normalize  R  to  unity  in  equation  (4),  the  Green’s 
function  equation  reduces  to 
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V;G(r  |  r0)  +  R2y2(r)G(r  |  r0)  =  S(r-r0).  (5) 

If  we  use  iterative  perturbation  theory,  the  zeroth  order  approximation  G(ol(r  |  ro)  is  obtained  by 
substituting  y  =  0 ,  and  is  governed  by  the  differential  equation 

V^G(0) (r  |  r0)  *  ^(r-r0).  (6) 

By  iteration,  we  get  a  differential  equation  governing  the  first-order  approximation  GU)(r  |ro), 
which  is  expressed  as 


V2G(1)(r  |  r0)  *  S(r-ro)-R2y2(r)Gm(r  |  r0).  (7) 

The  solution  to  equation  (6)  is  the  well  known  Green's  function  to  Laplace's  equation  in  polar 
coordinates  and  is  given  by  [6] 


G  (r  lro)  =  T— ln| 

4  n 


P2  +Po  -2 PPo  cos(0-0o ) 


\  +  P2Po  -IpPo  cos(£-£0) 


(8) 


Thus,  G'"(r  |  ro)  can  be  written  as 

G(1,(r  |  r0)  =  £  JS'G'V  |  r0)[S(r  -ro)-i?V2(r’)G(0,(r'  |  ro)]. 
Equation  (9)  can  be  simplified  to 

G(1)(r  |  O  =  G(0)(r  |  r0) -£,^'/?V2(r')G(0)(r'  |  r)G(0,(r'  |  r0). 


(9) 


(10) 


We  assume  that  y2( r)  is  piecewise  constant  in  6 ,  with  possible  discontinuities  at  ^=0,  nl 2,  n 
and  3tt/2  .  We  also  set  ro  =  0  (po  =  O,0o  =  0) .  As  a  result,  G(1>(r)  becomes  [1] 


G(1)(r)  =  —In  p- 
2  n 


R2 

8k  ' 


p1  +  p  -  2 pp  cos(6  -  &  ) 
1  +  p2 p  -2 pp  cos {0-6) 


p\np.  (11) 


Above,  y  =  y(  v)  for  r  in  the  ^th  quadrant.  Interchanging  the  order  of  summation  and  integration 
operations,  equation  (11)  leads  to 


R2 

8  ft2 


GU)(  r)  =  —In  p- 

2ft 


El/pj 


qn  /  2 
{q-\)7r  /  2 


dO  y2  In  p  In 


P  +P 


1 +  P2P 


-2 pp  cos (0-0  ) 
—  2 pp  cos {0-0  ) 


(12) 
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GU)(r)  is  an  approximate  expression  for  the  volume  Green's  function  for  Helmholtz  equation, 
including  terms  up  to  the  first  order.  Using  Green’s  theorem  [6]  and  imposing  Dirichlet  boundary 
conditions,  one  can  show  that  a  particular  component,  say  Ax  of  the  vector-potential  given  by 
equation  (3)  at  any  point  ro  is 

Ax(r0)  =  £[m/c.VrG(1)(r  |  rJK(r)  +  £</SG(1)(r  |  r0) IpJ .  (13) 

The  first  integral  in  equation  (13)  is  over  the  circumference  of  a  unit  circle  centered  atro.  The 
second  integral  in  equation  (13)  is  over  the  surface  of  the  same  unit  circle.  Now, 
hdc-Vr  — >  pddd  /  dp  at  p=  1  in  polar  coordinates.  Taking  the  derivative  of  G<u(r)  yields 


cG 


(i) 


dp 


'  p= i 


--4zf vr2  de'  . 

In  4n~^Jo  hq-\).n  p  —  2yC>  cos($  —  $  )  + 1 


(14) 


The  vector  quantity  VrGU)(r  |ro)  is  the  "surface  Green's  function",  and  equation  (13)  together 
with  equation  (14)  form  the  basis  of  a  floating  random-walk  algorithm  for  solving  a  scalar 
Helmholtz  equation  subject  to  Dirichlet  boundary  conditions.  The  random-walks  start  at  a  point, 
where  we  need  to  estimate  a  particular  component  of  the  vector-potential  given  by  equation  (3). 
The  random-walks  propagate  as  “hops”  of  different  sizes  from  circle  centers  to  circumferences, 
consistent  with  a  statistical  interpretation  [2]  of  (13).  Maximally  sized  circles,  subject  to 
limitations  imposed  by  iterative  perturbation  theory,  are  used  with  hop-location  probability  rules 
again  consistent  with  (13).  We  define,  with  each  hop,  a  numerical  weight  factor  derived  from 
(12).  The  product  of  these  weight  factors  over  a  walk,  multiplied  by  the  solution  at  the  problem 
boundary — where  the  walk  must  terminate — gives  a  statistical  estimate  for  the  component  of  the 
vector-potential  at  the  random-walk  starting  point.  All  this,  again,  is  entirely  consistent  with  a 
statistical  interpretation  of  (13).  We  can  thus  obtain  an  accurate  statistical  estimate  for  the 
component  of  the  vector-potential  by  averaging  over  a  statistically  large  number  of  random- 
walks.  Mathematically,  we  can  write  such  an  estimate 


A 


1  N 

—  TA 
N  n= 1 


(«) 


(15) 


where  N  is  the  number  of  walks  and  Au,)  is  the  nth-walk  estimate.  At  this  point,  it  is  easy  to  see 

the  relative  advantages  of  this  algorithm.  There  was  no  need  of  any  discretization  of  the  problem 
domain  leading  to  lower  memory  requirements  in  comparison  with  discretization  based 
algorithms.  It  can  also  be  noted  that  the  random-walks  are  independent  of  each  other,  which 
makes  the  algorithm  completely  parallelizable.  In  the  next  section,  we  will  describe  the 
benchmark  problems  that  have  been  solved  in  this  work. 
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4.  Homogeneous  Benchmark  Problems 


Two  materially  homogeneous  problems  have  been  chosen  in  an  insulating  and  conducting 
material  domain  respectively  for  a  problem  geometry  shown  in  figure  1.  In  each  of  these 
problems,  the  values  of  the  x-component  of  the  vector-potential,  Ax  at  gray-region  boundaries 

are  zero  at  top,  bottom,  left  and  right  edges.  We  have  chosen  a  forcing  function  that  varies 
sinusoidally  in  the  x  and  y  direction.  The  motivation  behind  such  a  forcing  function  is  that  any 
arbitrary  forcing  term  (that  is  piece-wise  continuous)  can  be  decomposed  into  an  infinite  number 
of  sinusoids  with  the  help  of  a  Fourier  expansion.  The  analytical  and  random-walk  results  are 
plotted  for  solution  points  along  the  centerline  starting  with  x  =  0  at  the  left  edge  in  figure  2  and 
figure  3  and  they  are  seen  to  be  in  good  agreement.  Table  1  summarizes  our  computational 
results.  From  Table  1,  it  can  be  noted  that  our  predicted  statistical  error,  is  consistent  with  the 
mean  absolute  error. 


Figure  1:  2D,  TM  homogeneous  test  problem,  the  geometry.  The  computational  region  is 
represented  in  white,  and  the  boundary  regions  where  vector  potential  is  known,  are  represented 
in  gray.  A  =  0  in  the  top,  bottom,  left  and  right  boundary  regions.  A  sinusoidal  forcing  tenn  is 
applied  in  x  and  y  directions:  / (x,  v)  =  sin(4;zx  /  L)siv\ny  /(2 L)\. 
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Ax 


x 

Figure  2:  2D,  TM  homogeneous  problem  in  insulating  medium.  Square  domain  with  L  =  100  pm, 
sr  =  2.7.  The  solid  line  represents  the  exact  analytical  solution  at  co  =  100  GHz.  The  dots  represent 
random- walk  solution  points  at  co  =  100  GHz.  A  normalized  length  scale  has  been  used. 

RefAj 


x 


Figure  3  a:  2D,  TM  homogeneous  problem  in  conducting  domain,  the  real  part  of  the  solution. 
Square  domain  with  L  =  10  pm,  p  =  1.8  pQ-cm.  The  solid  line  represents  the  exact  analytical 
solution  at  co  =  1  GHz.  The  dots  represent  the  random-walk  solution  points  at  co  =  1  GHz.  A 
normalized  length  scale  has  been  used. 


April  19-23,  2004  -  Syracuse,  NY  ©  2004  ACES 


20th  Annual  Review  of  Progress  in  Applied  Computational  Electromagnetics 


Im  (Ax) 


x 

Figure  3b:  2D,  TM  homogeneous  problem  in  conducting  domain,  the  imaginary  part  of  the 
solution.  Square  domain  with  L  -  10  pm,  p  =  1.8  pQ-cm.  The  solid  line  represents  the  exact 
analytical  solution  at  a>  =  1  GFlz.  The  dots  represent  the  random-walk  solution  points  at  co  =  1 
GFlz.  A  normalized  length  scale  has  been  used. 


Problem  # 

/(GHz) 

Random 

Walks 

Per  Point 
(xlO3) 

Mean 

Statistical  error 

Mean 

Absolute  Error 

Homogeneous 

Dielectric 

400 

10 

0.02 

0.02 

Homogeneous 

Conductor 

1.0 

10 

0.02+0.01/ 

0.01+0.01/ 

Table  1 :  Numerical  results  for  the  TM  homogeneous  benchmark  problems. 


April  19-23,  2004  -  Syracuse,  NY  ©  2004  ACES 


20th  Annual  Review  of  Progress  in  Applied  Computational  Electromagnetics 


5.  Conclusion  and  direction  of  future  work 

Summarizing,  we  have  tested  our  iterative-perturbation-theory  based  2D  floating  random-walk 
algorithm  on  a  homogeneous  TM  problem.  Excellent  agreement  has  been  obtained  between  our 
numerical  results  and  known  analytical  solutions.  Extension  of  this  algorithm  to  a  2D 
heterogeneous  TM  problem  as  well  as  the  development  of  a  3D  version  of  this  algorithm  will 
form  the  basis  of  future  work.  We  believe  that  with  additional  development,  this  algorithm  can  be 
useful  in  the  electromagnetic  analysis  of  complex  IC  interconnect  structures. 
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